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ABSTRACT 

Strong gravitational lensing of distant supernovae (SNe), particularly Type la's, has some exploitable properties 
not available when other sorts of cosmologically distant sources are lensed. One such property is that the "standard 
candle" nature of SN at peak brightness allows a direct determination of the lensing magnification factor for each 
well observed image. Another is that the duration of a SN event is of the same order as the differential time 
delays between the various lens images for roughly galaxy mass lensing objects. A relatively precise constraint 
on each image's magnification leads to better constraints on the lens mass model than are available in more 
familiar lens systems, and the comparable time scales of the photometric event and the time delay invite a variety 
of applications, including high precision measurements of the delay and the targeting of especially interesting 
phases of the explosion (including its very early stages) for intensive observation when they appear in trailing 
images. 

As an initial exploration of these possibilities we present calculations of SN lensing statistics in a "concordance 
cosmology" assuming a simple spherical model for lens mass distributions. We emphasize magnification and time 
delay effects. Plausible SN surveys, such as the proposed SNAP space mission, would discover several to some 
tens of strongly lensed SNe la per year, and at least a few of these will be at redshifts well beyond those that would 
be accessible via unlensed events. The total number of such anomalously high redshift SNe la will be a useful test 
of high redshift star formation models. SN surveys of finite duration will, of course, miss the appearance of some 
images, and the effect becomes large when the delay approaches the survey duration; we quantify this selection 
bias. Finally, we investigate how well the appearance of trailing images can be predicted based on various amounts 
of available information on the lensing event. Knowledge of the magnification factor for the leading (and brighter) 
image makes it possible to predict the appearance of a trailing image relatively accurately if the lens redshift is 
also known. 

Subject headings: cosmology: theory — gravitational lensing — supernovae: general 



I 1. INTRODUCTION 

Recent systematic surveys of distant type la supernovae 
(hereafter SNe la) strongly suggest the presence of dark energy 
which may dominate the total energy density of the universe 
(Riess et al. 1998; Perlmutter et al. 1999). The reason that this 
1 rather surprising conclusion is taken so seriously stems from the 
' fact that the SNe la are excellent (albeit empirical) "standard 
, candle" distance indicators, after an appropriate correction for 
the peak luminosity dependence on the shape of the individual 
lightcurve (Phillips 1993; Riess, Press, & Kirshner 1996). The 
exploitation of SNe la as cosmological probes has already been 
extensive (Riess et al. 1998; Perlmutter et al. 1999), and plans 
for yet more extensive observational studies, such as LSST^ 
(Large-aperture Synoptic Survey Telescope), are being rapidly 
developed. The most ambitious of these is the proposed satellite 
SNAP^ (SuperNova/ Acceleration Probe) which would gather 
~ 2000 SNe la per year by frequently imaging ^ 20 square 
degrees of sky. 

In this paper we explore some other ways in which the unique 
transient "standard candle" properties of SNe la might be ex- 
ploited for those few which happen to be strongly gravitation- 



ally lensed by an intervening object of roughly galactic mass. 
For the sake of specificity we will adopt the observational pa- 
rameters associated with the proposed SNAP mission unless 
otherwise stated. 

Very roughly 0. 1 percent of sources at z > 1 are expected to 
have multiple images due to strong gravitational lensing (e.g., 
Turner, Ostriker, & Gott 1984), and SNAP therefore would find 
at least ^ 2 lensed SNe la per year (e.g., Holz 2001). They 
will be qualitatively different from the lensed systems so far de- 
tected (e.g., ^ 60 QSO multiple-image systems^). First of all, 
the time-delay between the multiple images could be robustly 
and accurately determined for each object, at least in principle. 
Indeed, for a typical cosmological lensing time delay, of an or- 
der of a year, it is unlikely that one will observe all of the mul- 
tiple images simultaneously due to the finite duration of SNe 
la a month). This also means that some lensing events may 
be missed because the observation time is finite; one or more 
multiple images may appear only before or after the observing 
season or program. This is in marked contrast to QSO multiple- 
image systems where the images are observed simultaneously 
and their presence and geometrical arrangement in the sky is of- 
ten the chief indication of lensing. Furthermore, since the SNe 
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la are believed to be a reliable "standard candle", the magni- 
fication factor of the lensing can be determined directly. This 
is not feasible for any other sources because an object's intrin- 
sic luminosity is basically unknown. This can provide valuable 
additional information on the lensing potential. With such in- 
formation, one may indeed predict the location and the epoch 
of the additional trailing images, at least statistically. 

In the present paper, we first analytically calculate the ex- 
pected number of lensed SNe la to be found by a SNAP-liks sur- 
vey, including the effects of time delay bias due to a finite ob- 
servation/survey duration. Since lens systems with wider sepa- 
rations have larger time delays on average, time delay bias be- 
comes more significant for bigger image separations. Next we 
describe ways to predict the appearance of additional images 
of lensed SNe la. We then briefly consider the consequences 
of using a different and perhaps more realistic model for lens 
mass distributions. In conclusion, we discuss some of the im- 
plications and possible extensions of this work. 

Since we are not here concerned with determining cosmo- 
logical parameters via lensing (although this constitutes one of 
the primary purposes of the SNe la survey), unless otherwise 
specified we simply adopt the lambda-dominated universe with 
(fJo, Aoj/i) = (0.3,0.7,0.7), where flo is the density parameter, 
Ao is the cosmological constant, h is the Hubble constant in 
units of 100kms"^Mpc"\ the so called concordance cosmology 
(Ostriker & Steinhardt 1995; BahcaU et al. 1999). 



extinction (Steidel et al. 1999). The third model (SF3) is con- 
sidered because it has been suggested that the rates at high z 
may have been severely underestimated due to an unexpectedly 
large amount of dust extinction (e.g., Blain et al. 1999). These 
star formation rates are easily converted to those in different 
cosmologies (Porciani & Madau 2001), and are considered to 
span the range of reasonably realistic possibilities. 

The number rate of SNe la which occur between z and z+dz 
is then 

dN 7?SNeIa(z) ^ADliz) 



:(l+zr 



(5) 



dz l+z H(z)(l + z) 

where is the solid angle of the observed region and D^iz) is 
the angular diameter distance. 

3. LENSING STATISTICS 

3.1. Image Separation and Time Delay Probability 

Distribution 

In this discussion, we omit the magnification bias (Turner 
1980; Turner et al. 1984) because the intrinsic luminosity func- 
tion at peak brightness of Type la SNe is quite narrow^. We 
consider that magnification of SNe at redshifts beyond those 
normally accessible to the survey, another form of magnifica- 
tion bias, separately in the next subsection. 

We now consider lensing objects with a the Singular Isother- 
mal Sphere (SIS) density profile: 



2. NUMBER COUNTS OF SUPERNOVA lA 

The number of SNe is closely related to the star formation 
rate /?sf in the universe. For a Salpeter initial mass function 
(^imf(M) 3c O.1M0 <M< 125Mq), the rate of Type la 

events /?sNeia is calculated as follows (Madau, Delia VaUe, & 
Panagia 1998) 

^SNela(z) = 

vlt^RsF{z{t'))dt'jl'^'' exp [-(t-t'-tM) /t] 4>mp{M)dM 

r / M(j)mi^{M)dM 

where Mc = max[3M0,(lOGyr/f')°'*M0] is the mmimum mass 
of a star that reaches the white dwarf phase at time f ' (assuming 
all systems with the primary star of mass 3M0 <M < SMq are 
possible progenitors of SNe la), tu = lOGyr/(M[M0])^^ is the 
standard lifetime of a star of mass M, r is a characteristic ex- 
plosion time scale (corresponding to an effective "delay" time 
scale after the collapse of the primary star to a white dwarf), rj 
is the SNe la explosion efficiency. We adopt three representa- 
tive star formation rates used by Porciani & Madau (2001). The 
explicit forms of these in an Einstein-de Sitter universe are 

exp(3.4z) , , 

^SFi(z) = 0.462/1 MQyr-'Mpc-^ 



exp(3.8z) + 45 

RsFiiz) = 0.230h^^j^^MQyr-'Mpc-\ 
exp(3.4z) + 22 

exp(3.05z-0.4) , 

Rspiiz) = 0.308/i . ,, M0yr"'Mpc" 



(2) 



(3) 



(4) 



exp(2.93z)+15 

The first model (SFl) includes a correction for dust redden- 
ing, and matches most measured UV-continuum and Ha lumi- 
nosity densities (Madau & Pozzetti 2000). The second model 
(SF2) is also possible because of the uncertainties associated 
with the incompleteness of data sets and the amount of dust 



(6) 



2-iTGr^ ' 

where v is a one-dimensional velocity dispersion and define the 
characteristic scale length 

2 Z)oL^?LS 



Co = 47r 



Dos 



(7) 



where Z)ol> Dqs, and Dls are the angular diameter distances 
to the lens, to the source, and between the lens and source, re- 
spectively. Then the lens equation becomes y =x-x/|x|, where 
X and y are image and source positions in each plane normal- 
ized by ^0 and (.oDos/Dqi^, respectively. This has two solu- 
tions x±=y±l, and each image will be magnified by a factor 
/i± = (1/y) ± 1. The angular separation of two images is given 
by 

6(x.-x-)^_g /vyA^_ (8) 
\C/ Dna 



0-- 



DoL vc/ Dos 

The differential time delay between two images can also be cal 
culated as 

/DolDls 



cAf (y) = 



: 327r" ( - 



£>os 



-(l+ZL)y. 



(9) 



The cumulative and differential probability distributions of 
strong lensing are 



P{>6;zs) 



= dzj 

JO Jv„, 



' cdt , 
c^sisi— (1 + zl) <A(v), 



p{e-zs) = --7^p{>e-zs). 

da 



(10) 



(11) 



where crsis = tt^^ is the cross section of strong lensing, v^in = 
{9 /%T^fl^{Dos/Di^syi^c, and 4>{v) is the velocity function of 
lens galaxies. We adopt the velocity function: 



(j)(v)dv -- 



1±. 
InlO 



exp 



'V2v 



fi/2S 



dv 

— , (12) 

V* 



^ This is equivalent to assuming that the siu-vey monitors every field sufficiently frequently to catch each SN event near its peak brightness, as is the case for SNAP. 
The bias could not be neglected for a more traditional survey which finds most objects well after their maximum brightoess. 
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Fig. 1 . — Supernova rates /?sN (eq- [1]). number rates of SNe la dN /dz (eq. [5]), and number rates of lensed SNIa dN\cns/dz = / dN / dzP(z)dz as a function of 
redshift of SNe la, on the basis of the model described in §2. In the middle and bottom panels, observations of 20 square degrees field is assumed and no magnitude 
limit is imposed. For RsN> four observed SN la rates are also plotted by filled triangles; z ~ 0.01 (Cappellaro, Evans, & Turatto 1999), z ~ 0. 1 (Hardin et al. 2000), 
z ~ 0.4 (Pain et al. 1996), and z ~ 0.55 (Pain et al. 2002). The efficiency parameter r) is adjusted so as to reproduce the observed SN la rate at z ~ 0.01. The 
"concordance" lambda-dominated universe with (Q,Q,XQ,h) = (0.3,0.7,0.7) is assumed. 



where = 7.3 x 10"2/^^Mpc-^ (3 = -1.3, and = 247kms-'. 
This distribution function is based on the Southern Sky Redshift 
Survey and is derived by Gonzalez et al. (2000). Averaging the 
above probability function for fixed zs over the observed rate of 
SNe (eq. [5]) yields the number rate of strong lensing; 



dN 

PiO;zs)—dzs. 
azs 



(13) 



Similarly, the joint probability distributions of time delays 
and image separations are (see eqs. [29]-[31] of Oguri et al. 
2002) 



P{>At,e;zs) 

dZL 



dv 
M 



cdt 



d(At) 



dZL 



P(> At,9;zs), 



C14) 



(15) 



where N^(> At) = 1 -y^j^ and jmin can be calculated from 
Af = Af(yniin)- These joint probability distributions divided 
by P{6;zs) give the conditional probability distributions P(> 
At\0-Z5)and P(At\0-zs). 

Figure 1 shows the predicted number rates of SNe la and also 
the expected number rate of lensed SNe la. The number rates 
of SNe la are calculated from equation (5). It is clear from 
this figure that the number of SNe la strongly depends on the 



star formation rate and its evolution, as indicated by Madau et 
al. (1998). The number of lensed SNe la also shows large dif- 
ferences between models, reflecting the above sensitive model- 
dependence of the number of SNe la. Our calculation predicts 
that the number of strongly lensed SNe la will be between a 
few and a few tens per year for a SNAP-liks survey. These lens- 
ing rates are roughly consistent with those calculated by Wang 
(2000) and Holz (2001). 

3.2. SNe beyond the Magnitude Limit 

Gravitational lensing magnifies SNe, thus some SNe which 
exceed the magnitude limit if they are unlensed might be ob- 
served due to their magnification (e.g., Kolatt & Bartelmann 
1998; Porciani & Madau 2000; Sullivan et al. 2000; Goobar et 
al. 2002). If we neglect the effect of gravitational lensing, the 
apparent magnitude of SNe at peak can be expressed as 

mx =MB + 51og(DL(z)[Mpc]) + 25 + A:Bx, (16) 
where mg = -19.4 is the peak magnitude of SNe la, Z)l(z) = 
(1+z)^Da(z) is luminosity distance, and Kbx is the single- 
or cross-filter K-corrections. Therefore, even for high-z SNe 
la whose unlensed apparent magnitude exceeds the magnitude 
limit, nix > mum, they are observed if magnified by a factor /i 
satisfying 

/i>10°'4<'«''-'"-^=Ai,nin. (17) 

We consider the following two cases: (1) /i_ > /imin- This cor- 
responds to the case that both lensed images are observed. (2) 
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/x+ > /Xniin- This means that at least one image is observed. In 
each case, the cross section for lensing is 

aifi^> ^min) = 7 — ■ (18) 

We calculate the numbers of such lensed SNe la, and the re- 
sults are shown in Figure 2. We assume the magnitude limit is 
wiim = 30, and also impose a requirement that the photometry 
must extend to 3.8 magnitudes below peak. We take account 
of the intrinsic dispersion of SNe la peak magnitudes assuming 
the Gaussian distribution with dispersion a,„ = 0.15 (Porciani 
& Madau 2000). As seen in the figure, the expected number of 
such lensed SNe la depends very strongly on models of star for- 
mation history. In most models, however, more than one high-z 
(z ~ 3) SNe la per year is expected to be observed. Actually 
the difference simply comes from the difference in supernova 
rates, thus we can infer supernova rates at high redshifts from 
the observed number of lensed SNe la. The supernova rate as 
a function of redshift is useful not only to constrain progeni- 
tor models and star formation history (e.g., Yungelson & Livio 
2000) as shown in these plots, but also to test other possible 
redshift dependence of the SNe la rate (Kobayashi et al. 1998). 
The SNe la rate also depends on the assumed cosmological 
parameters. The cosmological parameters flo and Aq are now 
determined with ~ 10% accuracy in the combined analysis of 
Cosmic Microwave Background (CMB) experiments, SNe la, 
and large scale structure surveys (e.g., Sievers et al. 2002), and 
we assume the ilo uncertainty to be ilo = 0.3 ± 0.05 in the flat 
universe (17o + Ao = 1). We also plot the resulting numbers of 
lensed SNe la for the constant /?sn model due to the uncer- 
tainty. We find that this level of uncertainty in ilo does not sig- 
nificantly change the expected number of SNe so that we could 
still distinguish between models of the star formation history. 
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Fig. 2. — Numbers of lensed SNe la whose unlensed appai'ent magni- 
tudes are fainter than the magnitude limit niiip, = 30 for the observed region 
of f! A = 20 deg^ . Lines ai'e same as in Figure 1 . For the constant Rsn model, 
we show the effect of different cosmological parameters (f2o = 0.3 it 0.05, 
^lo + X() = 1) by the same three lines and shadings. We also impose a require- 
ment that the photometry must extend 3.8 magnitudes below peak brightness 
(thus effectively mij,,, = 26.2), and K-corrections are neglected. The limiting 
peak magnitude ninn, = 26.2 roughly corresponds to z ~ 1.7. The scatter of 
peak magnitude, am = 0.15, is taken into account. A lambda-dominated uni- 
verse is again assumed. 



If both images are observed, we can infer the intrinsic mag 
nitude of SNe la. For example, from the ratio of luminosities, 

r= — 

the magnification factors become 



(19) 



2r 
7^' 



r-l 



(20) 
(21) 



The magnification factors also can be derived from the im- 
age separation and differential time delay if the redshift of the 
lens object is measured (see eq. [30]). The reconstruction of 
magnification factors could be used, in principle, to derive the 
distance-redshift relation at high-z and to estimate the cosmo- 
logical parameters more robustly. 

4. TIME DELAY BIAS 

Since there is a time delay between multiple images, strong 
lensing statistics of transient phenomena such as SNe inevitably 
involve some missing events due to the finite duration of the sur- 
vey observations. Therefore, in strong lensing statistics of SNe 
we should take account of this "time delay bias", especially for 
large image separations. 

For the joint probability distributions, the time delay bias is 
included as follows 



P''^(Af,0;zs) = P(Af,0;zs)/(Af), 



(22) 



P™(> Af,0;2s)= / PiAt',9;zs)fiAt')d(At'), (23) 

J At 

where /(Af ) is the fraction of lenses with time delays At that 
can be observed (the superscript TB explicitly indicates the 
Time-delay Bias). The image separation distribution then be- 



P^''(0;zs)= / P(Af',0;zs)/(AfV(Af'). 



(24) 



Then the corresponding conditional probability distributions 
are computed as 

P™(At,e;zs) 



P'^(At\e;zs) = ' 



(25) 



P™(0;zs) ' 

using P™(0;zs) instead of P{6;zs) which does not take account 
of the time-delay bias (eq.[ll]). 

If the observational monitoring is carried out continuously 
for a period of fobs, for instance, f{At) is given by 

'' Af , 

1- — (Af <fobs), 
*obs 

(Af > fobs). 



/(Af) : 



(26) 



Then P™(6';zs) reduces to 

r 1 i-'ok 
P™(0;zs)= 1 / Pi> At\9;zs)d(At) 

fobs Jo 

= T(d,zs)P(0;zs). 



P(0;zs) 



(27) 

This means that the time delay bias for the image separation 
probability distribution is simply expressed by the multiplica- 
tion factor T(9,zs)- 

Probability distributions of strong gravitational lensing in- 
cluding time delay bias (eq. [27]) are shown in Figure 3, 
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Fig . 3 . — Probability di.stributions of strong gravitational lensing. The redshift of SNe is fixed to zs = 1 or 2 (right). The time delay bias (§4) is included, and 
continuous observations of t^bs = c», 3yr. and lyr are considered. These probabilities are plotted for the three cosmological models with h = 0.7; (Qq.Xq) = (0.3,0.7), 
(0.3,0.0), and (1.0,0.0), in solid, dashed and dotted lines, respectively. 



where we consider three cosmological models with h = 0.7; 
(Oo, Ao) = (0.3,0.7), (0.3,0.0), and (1.0,0.0). These probability 
distributions have been used to constrain cosmological constant 
(Turner 1990; Fukugita et al. 1992; Kochanek 1996; Chiba & 
Yoshii 1999), and are indeed useful as an independent mea- 
surement of cosmological parameters in the case of SNe sur- 
vey (Wang 2000; Holz 2001; Goobar et al. 2002). These plots 
show that the time delay bias is more important for larger 6, 
because the time delay At is larger on average as 6 increases. 
The amount of the time delay bias is, however, quite small for 
lensing of typical an image separation of 9 ^ 1", if the obser- 
vation time is larger than 1 year. For larger separation lensing, 
9 ~ 3", the time delay bias changes lensing probabilities by a 
factor 2 and thus is important in quantitative discussions. 

Since the time delay bias is more important for larger sep- 
aration lensing, we also calculate probability distributions for 
wide separation lensing in the lambda-dominated Cold Dark 
Matter (CDM) model, assuming a fluctuation amplitude cts = 1 
for simplicity. Although the expected wide separation lensing 
rate due to CDM halos is much smaller than for galaxy lensing, 
we still expect that wide separation lensing will be observed be- 
cause core-collapse SNe as well as SNe la can be used for lens- 
ing statistics, which greatly increases the number of SNe (e.g., 
Goobar et al. 2002). Wide separation lensing is expected to re- 
flect the properties of dark halos rather than (the visible parts 



of) galaxies, thus it has been used to constrain the abundance 
of dark halos (Kochanek 1995) and density profile of dark ha- 
los (Maoz et al. 1997; Wyithe, Turner, & Spergel 2001; Keeton 
& Madau 2001; Takahashi & Chiba 2001; Li & Ostiiker 2002; 
Oguri et al. 2002). We adopt the generalized form (Zhao 1996; 
Jing & Suto 2000) of the density profile proposed by Navarro, 
Frenk, & White (1997, hereafter NEW): 

Peril 



p{r): 



(28) 



(rA,)"(l+r/r,)^-"' 

where = ryii /cvii and Cyir is the concentration parameter We 
adopt the mass and redshift dependence reported by Bullock et 
al. (2001) for a = 1, and generalize it to a ^ 1 by the multiplica- 
tive factor (2- a) (Keeton & Madau 2001). We also take ac- 
count of scatter of the concentration parameter which has a log- 
normal distribution with a dispersion of CTc = 0.18 (Jing 2000; 
Bullock et al. 2001). The characteristic density 6^ can be com- 
puted using the spherical collapse model (see Oguri, Taruya, & 
Suto 2001). As the mass function of dark halos, we adopt the 
fitting form derived by Sheth & Tormen (1999). From these, 
we predict probability distributions of wide separation lensing 
taking account of the time delay bias (see Oguri et al. 2002 
for the calculation of time delay probability distributions in the 
case of generalized NEW density profile), and results are shown 
in Figure 4. Here we focus on large separation lensing (6 > 5") 
because the relation between small and large separation lensing 
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0[arcsec] 

Fig. 4. — Probability distributions of strong lensing at wide separation, calculated by assuming the generalized NFW density profile (eq. [28]) and Sheth & 
Tormen (1999) mass function. We adopt the cold dark matter model in a lambda-dominated universe with the mass fluctuation amplitude erg = 1. The time delay 
bias (§4) is included. 



depends strongly on the model of galaxy formation and is dif- 
ficult to determine unambiguously (Oguri 2002). This plot in- 
dicates that time delay bias is quite significant; it can suppress 
the lensing probability by one or two orders of magnitude. It is 
also found that the suppression due to time delay bias is larger 
for a = 1 .5 than a = 1 .0 because for a fixed separation 9, the 
time delays in the case of a = 1 .5 are on average larger than 
those in the case of a = 1.0 (Oguri et al. 2002). This slightly 
compensates for the difference of lensing probabilities between 
various values of a, but the difference is still large (about one 
order of magnitude between a = 1.5 and a = 1.0). Therefore 
we conclude that statistics of wide separation SNe lensing still 
could provide a useful probe of density profile. 

We also show the effect of the time delay bias on the time 
delay probability distribution in Figure 5. As seen in the figure, 
the time delay probability distributions for large 9 are strongly 
affected by the finite duration of observations. 

5. PREDICTIONS FOR TIME-DELAYED (TRAILING) IMAGES 

Since the primary purpose of the proposed SNe la surveys 
is to construct the Hubble diagram, one may always assume 
that the source redshift of each SNe la, Zs, is determined spec- 
troscopically. We are interested in unusually bright SNe la in 
the diagram. If the lensing object is approximated by an SIS, 
there are two images, x±, and the brighter image (x+) arrives 
first. Thus the difference between the observed magnitude of 



those outliers (with respect the typical magnitude of the SNe 
la) and the average magnitude at the redshift can be ascribed to 
the lensing magnification fi+. While the standard deviation of 
the SNe la peak magnitude corrected for the lightcurve shape 
method is typically O.lSmag (Porciani & Madau 2000), multi- 
ple lensing images are produced if > 2 or equivalently more 
than 0.75mag. Thus Scr outliers are strong candidates for mul- 
tiple image lensed SNe la. 

The lensing object (most likely an early-type galaxy) for 
those candidates is typically half way out in affine or angular 
diameter distance. For instance, an L* galaxy at z = 0.5 is about 
24th magnitude in the B-band, and the surface number density 
at this magnitude limit is ^ 30 per arcmin^ or so. Since the typ- 
ical image separation 9+ is 1", the average number of galaxies 
within that angular separation from each lensed SNe la is much 
less than 1 . Therefore it should be relatively easy to find can- 
didates for the lensing galaxy and thus determine 6+, and even 
possibly the redshift of the candidate galaxy zl- 

Our next goal is to predict the location of the second image 
9-, and the time delay At by which the second image trails the 
first and brighter one. We consider the following three cases: 
(1) the magnification factor the redshift of the lens object zl, 
and the angular separation from the lens object 9+ = £,ox+/Do]^ 
are all known; (2) only /i+ and 6+ are known; (3) just ^+ is 
known. In doing so, we have to assume a specific model for the 
lensing potential. We mainly show results for the SIS lensing 
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At[yr] 

Fig. 5. — Conditional probability distributions of differential time delays. The time delay bias (§4) is included (as in Fig. 3). A lambda-doininated universe is 
assumed. 
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Fig. 6. — Distributions of lens redshift zl (eqs. [32] and [38]). Thick hnes are probabiUty distributions of zl when 8 is unrestricted, P{zl) (eq. [32]) while thin 
lines are those when 8 is fixed, P(zh\S) (eq. [38]). A lambda-dominated universe is assumed. 



model (eq. [6]) but also present a comparison with the NEW 
halo model (eq. [28]). 

5.1. Case 1: fi+, 9+, and zl are known 

We can obtain these three quantities when the lens candidate 
is identified and the redshift of that lens candidate is known. 
In this case, both the separation 6 and the time delay At are 
uniquely determined as 



^ = 29+'- 



and 



Af=l^^(l+ZL)02 1 



2c Dls 



1 



(29) 



(30) 



We now consider errors induced by observable quantities. 
We assume that the redshifts are measured spectroscopically, 
thus errors from redshifts are negligible. As for the image sep- 
aration 6+, SNAP would have an angular resolution of 0.1", 
which is also sufficient to determine 9+ accurately for most pur- 
poses. The most important source of uncertainty is therefore 
because the magnification estimates may be inaccurate due 



to substructure in the lens galaxies (Mao & Schneider 1998), 
dust extinction and/or the intrinsic spread in corrected SNe la 
peak luminosities, photometry, K-corrections and so forth. If 
11+ is sufficiently larger than 1 (this is correct for most strong 
lensing cases), we obtain Af cx (//+)"' from equations (29) and 
(30). This means that errors in /i+ directly affect Af ; e.g., a 20% 
error in /i+ results in a 20% uncertainty in the Af estimation. 

5.2. Case 2: /i+ and 9+ are known 

This is the case in which there is a lens candidate but the 
redshift is not known (yet). We can predict the image separa- 
tion from equation (29), but the time delay is ambiguous due 
to Zl, and has some probability distribution which reflects the 
distribution of Zh- We rewrite equation (30) as 
Af(^+-1) 1 Dol£>os, 



-{\+Zx} = F. 



(31) 



9^ 2c Dls 

Therefore we can obtain the probability distribution for the 
combination of Af(/x+- 1), instead of Af alone. To derive this, 
we should calculate the conditional probability distribution of 
Zl for fixed 9: 

Piz^\e;z.) = -^^^, (32) 
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Fig. 7. — Predictions for the trailing image in "Case 2" (§5.2). Probability distributions of (p+ — I) At (eq. [34] times (p+ — l)At) are plotted for various zs ™d 9- 
A lambda-dominated universe is assumed. 



where 



P(zlAzs)-- 



dvd(j)cdt 3 
-T^-p-}— (1+Zl) ctsis 
ai) av flL 



(33) 



We plot P{z]_\9;zs,) for various 9 values IN Figure 6. Then the 
distribution of Af(/i+- 1) is 



P(Af(M+-l)) = P(zL|0;zs) 



dF 
dzh 



1 



(34) 



The results are shown in Figure 7. These results are of course 
affected by observational uncertainties of fj,+ as mentioned in 
§5.1. If /i+ is sufficiently large, errors in do not affect 6 
since 9 ~ 29+ (see eq. [29]). On the other hand, Af is directly 
affected by ii+ because we obtain the probability distribution for 
A/(/i+- 1). The width of this probability distribution is, how- 
ever, about one order of magnitude, while errors induced by /i+ 
are probably only a factor of 2 or so (Metcalf & Madau 2001; 
Chiba 2002; Dalai & Kochanek 2002). Therefore errors in /i+ 
are not so serious as in the previous case. 



In practice, we find a fitting formula of f (Af (/i+- 1)): 
P(At(fi+-l))d(At(p+-l)) = 

,2' 



1 



Incr 



exp 



{ln(Af(/i+-l))-lnflf 



2^2 



fl = 4.28(-110+158z^ 



°-2i2)^i«(3.68 + 0)-°"52 



t/ln(Af(^+-l))(35) 
[day], (36) 



cr = 0.693 - 0. 1 15(ln6i) -I- 0.02 11 (In 6*)^ 

+0.00347(ln6i)^^ - 0.00106(ln6')^ (37) 
where 9 is in units of arcsec. This formula is valid for 0.5 < 
Zs < 4.0 and 0.1" < 6* < 10" which cover the range of our typ- 
ical interest. The accuracy is < 10% around the peak (within 

- 1.5cr). 

5.3. Case 3: fi+ is known 

This is the case in which a candidate for the lensing object is 
not identified. The distribution of zl for an unrestricted 6 is 

P(zl;zs)=-^ f dv^^d+ZLfasis, (38) 
P(zs) Jo dv dzL 
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where we normalize by the total probability of strong gravita- 
tional lensing: 



P(Z5)= / dZLP(ZL;Z5). 



(39) 



From this distribution of zs, we can calculate the probability 
distribution for At(p+-l)/9^: 



(40) 



Therefore we can predict the probability distribution for a com- 
bination of the time delay, the magnification and the image sep- 
aration. Figure 8 plots this probability distribution for various 
Zs- Errors in /i+ are not important for the same reason described 
in §5.2. 

Again we find a fitting formula for P(Af(/i+- l)/9^): 
P(At(^i^-l)/9^)d(At(p^-l)/9^) 



lira 



exp 



{ln{At(^i+-l)/9^)-lna} 



2a2 



xdln{At(fi+-l)/9^) , 



fl = -110+158z'^-2'2[day], 
(7 = 0.846. 



(41) 
(42) 
(43) 



This formula is valid for 0.5 ^ Zs 4-0, and the accuracy is 
< 10% around the peak (within ~ 1 .5a). 
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Fig. 8. — Predictions for the trailing image in "Case 3" (§5.3). Probability 
distributions of {fu-i-)At/9- (eq. [40] times {fu-l)At/e'^) are plotted for 
various zs ■ A lambda-dominated universe is assumed. 



5.4. Case 2 & NFW profile 

In this subsection, we consider a case in which the density 
profile of the lensing objects is well described by the general- 
ized NFW profile (eq. [28]) instead of by a SIS. We retain the 
velocity function of galaxies (eq. [12]) previously used, since 



we mainly want to elucidate the effect of different density pro- 
files. In this case, differential time delays and image separations 
are approximated as (Oguri et al. 2002) 

'^'"'''^^-(l+zDy, (44) 



Af : 



cDolDls 

£>0L 



(45) 



where Xi is a radius of the tangential critical curve normalized 
by ^0 = Ts. We also assume that /i+ is approximately given as 

(46) 



2 y 



where is the radius of the radial caustic. The explicit form of 
/ito may be found in Oguri et al. (2002). From these, 

^i+At = 9 (1+Zl), (47) 

4cxt Dls 

and thus we can derive the probability distribution of /i+Af. 
Figure 9 shows the probability distribution of /i+Af for inner 
slopes a = 0.5, 1.0 and 1.5 in dashed, dotted and thin solid 
curves, respectively. This figure indicates that a smaller a pre- 
dicts larger values of /i+Af. The dependence of Af alone on a 
has the opposite sign ; a steeper inner slope of the density have 
larger time delays on average (Oguri et al. 2002). The reason 
is that a shallower density profile tends to produce larger /i+ 
values and thus to cancel the dependence of Af and /i+ on a. 
The effect of varying the assumed lens density profile is thus 
significantly reduced by competing effects. 

6. DISCUSSION 

We have studied strong gravitational lensing of distant su- 
pernovae with particular attention to their magnification and 
time-delay statistics. Since supernovae are both "standard can- 
dles" at peak brightness and transient phenomena, unlike more 
conventional lensed sources such as quasars, they have some 
unique and relatively unfamiliar properties. 

One such feature of supernova lensing is that we can deter- 
mine the magnification factor directly from observations, not 
just the magnification ratios between different images. Since 
SNe la are known to be an excellent standard candle, unusually 
luminous supernovae are always strong candidates for strong 
lensing. We have shown that the location and the time de- 
lay of the trailing images of a lensed supernova with a given 
magnification factor can be predicted with useful accuracy. We 
have further considered several cases, depending on whether 
the lensing galaxy can be identified or not and whether its red- 
shift is known or not. Such predictions will allow targeted 
observing programs to study exceptionally interesting phases 
of SN explosions and the determination of extremely accurate 
time delays. While we have mainly considered lenses with an 
SIS density profile, we find that our results are not qualitatively 
sensitive to variations in this density profile, and specifically, 
that generalized NFW profile lenses produce quite similar ef- 
fects. 

Due to the finite duration of the event, strongly lensed SNe 
will not always have multiple images observable simultane- 
ously. In many cases, the second (usually fainter) image will 
appear after the first (brighter) image has faded away. This 
leads to an observational bias against the detection of multiple 
lensing images in any realistic SNe survey. We have calculated 
this time delay bias analytically, on the basis of time delay prob- 
ability distributions derived by Oguri et al. (2002). We find that 
time delay bias significantly changes the expected number of 
lensed SNe, especially at wide separations. More specifically, 
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if the observational survey lasts of an order of a year, the lensing 
probability is suppressed by more than one order of magnitude 
at ^ 10". The suppression of the lensing probability is greater 
for a steeper inner density profile in the lensing objects. This is 
simply because a steeper inner profile yields a larger time delay 
(Oguri et al. 2002). We note that the time delay bias may have 
less effect on ground based SN surveys which are expected to 
operate for a longer term, such as LSST or a SN pencil beam 
survey like that proposed by Wang (2000). 

In this paper we have studied spherically symmetric lenses. 
Although the inclusion of small ellipticities has little effect 
on the lensing cross section (Blandford & Kochanek 1987; 
Kochanek & Blandford 1987), it can yield lensing systems with 
four images. In such a more realistic case, our predictions for 
the time delays and image separations should correspond ap- 
proximately to those for the pair of images with the largest sep- 
aration in each quadruple (their fractional errors are expected to 
be of order the ellipticity). Time delays among the other images 
can be much smaller than our predictions; e.g., the time delay 
between A 1 - A2 in PG 1 1 1 5-1-080 is expected to be much smaller 
(of order of one hour) than the time delay between B-C (~ 25 
days) (e.g., Keeton & Kochanek 1997). The time delay bias for 
such systems might then be greatly reduced. In addition, four 



images usually appear when the position of the source is close 
to that of the center of the lens galaxy. This will also cause four- 
images systems to have systematically smaller time delays. We 
thus expect that the ratio of four image to two image lens sys- 
tems will be larger for lensed SNe than for other lens systems 
in which all images are continuously present. 

Another pleasing aspect of strongly lensed SNe la is that they 
allow one to make a fairly straightforward theoretical predic- 
tion for a cosmologically distant phenomenon that is then sub- 
ject to direct quantitative verification on a humanly practical 
time scale. This possibility is relatively rare in astronomy ex- 
cept in case involving intrinsically periodic phenomena, such 
as orbits or pulsar emission, where the "prediction" is basically 
empirical extrapolation rather than truly theoretical. It is partic- 
ularly rare in a cosmological context. Although one certainly 
would not expect major surprises in comparing such predictions 
to future observations, it is still an important opportunity to test 
and validate our basic understanding of cosmology and gravi- 
tational theory. 

This research was supported in part by the Grant-in-Aid 
for Scientific Research of JSPS (12640231, 14102004) and by 
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